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Abstract 

We consider integro-differential models describing the evolution of 
^ . a population structured by a quantitative trait. Individuals interact 

competitively, creating a strong selection pressure on the population. 
On the other hand, mutations are assumed to be small. Following the 
formalism of [17J, this creates concentration phenomena, typically con- 
sisting in a sum of Dirac masses slowly evolving in time. We propose 
a modification to those classical models that takes the effect of small 
populations into accounts and corrects some abnormal behaviours. 
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1 Introduction 

1.1 The first model 

We study the dynamics of a population subject to mutations and selection 
due to competition between individuals. Each individual in the population 
is characterized by a quantitative phenotypic trait x (for example the size of 
the individual, their age at maturity, or their rate of intake of nutrients). For 
simplicity, x is taken here in M even though all the arguments could easily 
be extended to higher dimensional cases. 
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Probabilistic models are usually considered as the most realistic in that 
setting. They consist in life and death processes for each individual Xj, 
Poisson processes more precisely, with for instance birth rate b{Xi) and a 
death rate which increases with the competition between individuals, for 
example 

When a birth occurs, it simply adds another individual with the same trait, 
except when a mutation takes place with small probability. In that case the 
new individual has a different random trait, obtained through some distri- 
bution K. We refer to [25] and the references therein for a nice introduction 
to the probabilistic approach. 

Of course this is only one possible model and there are many variants. One 
can modify the interaction between individuals for instance by introducing 
explicit resources (with chemostat like interactions maybe). The competition 
could influence both the birth rate and the mortality rate... 

When the total number of individuals is too large (it can easily reach 
10^° — 10^^ for some micro-organisms), it becomes prohibitive to compute 
numerically the solution to this process. In that case one expects to be 
able to derive a deterministic model as a limit of large populations. Such 
of derivation was proved in [8] and one obtains integro-differential equations 
hke 



dtu{t,x)= {b{x)-d{x)- / I{x-y)u{t,y)dy\ u{t,x) + M{u){t,x), (1.1) 

where the mutation kernel is for instance 

M{f){x) = f K{z) {h{x + z) f{x + z)- h{x) f{x)) dz. 



1.2 The scaling of fast reaction and small mutations 

Eq. (11.11) is easy and fast to solve. However in most situations, small param- 
eters appear and complicate that. This is a consequence of two small scales 
that are typical for this problem: 

• The rate of mutation is very small in comparison to the reproduction 
rate. As we wish to see the evolution of traits generated by the muta- 
tions, this means that we need to rescale the equation in time and will 
hence get a large reproduction rate. 



• The size of each mutation is sinalL 

To simphfy again our equations, let us now assume that 6=1 and denote 
r{x) = 1 — d{x) the reproduction rate of an individual without competition. 
Taking the two scalings into account, one has in fact to deal with an 
equation like 

dtUs{t,x) = -{r{x) ~ I-kUe) Ue{t,x) + Me{Ue){t,x), (1.2) 

where the mutation kernel now reads 

Me(/)(x) = - f K{z){h{x + ez)f{x + ez)-h{x)f{x))dz. 

Note that Ms{f) is indeed of order 1 if 6/ G C^ for instance. 

Eq. (11.21) is now much more delicate. The properties of the solution 
might depend on e (its smoothness for instance) and solving it numerically 
can become again very costly if e is too small (a typical value for many 
applications would be e ~ 10~^). Therefore one would wish to derive a new 
model as e — !■ 0. 

Eq. (II. 2p is strongly similar to a reaction-diffusion equation with a strong 
reaction term. However a crucial difference here is that the reaction term 
is non local. As we will see this completely changes the behaviour of the 
solution. 

This scaling was introduced in [17] and formally studied there. We briefly 
reproduce the main argument here. The starting point is to introduce a large 
deviation scaling (which makes sense in view of the original probabilistic 
interpretation of the model) 

iPeit,x)= elogUe. (1.3) 

Then it is easy to see that Eq. (11.21) becomes 

dt(fe = r- I-kUe + HsiiPe), (1-4) 

where 

H,{f) = f K{z) (e(/(-+^^)-/(-))/- - 1) dz. (1.5) 

Jr 

Using the theory of viscosity solutions to Hamilton- Jacobi equations, one can 
pass to the limit in (II. 4p and obtain for the limit (f of (fe 

dtip = r -Ii<u + H{d.^ip), (1.6) 



with u the weak hmit (in the space of measure) of u^ and 

H,{f) = [ K{z) (e(/(-+-^)-/(-))/- - 1) dz. (1.7) 



Of course at the hmit, if and u are no more connected by a relation hke (11.31) . 
Therefore Eq. (11.61) is no more closed and the question of how to recover u 
from ip is one of the main difficulty here, which is now fortunately better 
understood though. 



1.3 The problem with Eq. (O) 



Let us focus here on one other delicate issue with this approach which is the 
main motivation for the current work. In the scaling under consideration, 
one has growth or decay of order exp{C/e). In particular one can see that 
the ratio between the maximal value of Ug and the value at most other point 
is of this order exp(C/£:). However if we come back to the starting point, 
which means a total population of 10^*^ — 10^^ and e ~ 10~^, then there is 
an obvious problem. If (and it can be proved) u^ is of order exp{C/e) over a 
fixed interval of traits / then the total population over this interval includes 
in fact much less than one individual! 

This has several consequences. First of all it means that given the scales 
under consideration here, the limit of the probabilistic models to the deter- 
ministic one is not fully justified (the total number of individuals would have 
to be much higher with respect to e). This, in itself, could be ignored as any- 
way with 10^^ individuals the obvious solution, use the probabilistic model, 
is not practical. Hence one could reasonably accept to still work with the 
deterministic model if its predictions were qualitatively in agreement with 
the behaviour of the stochastic one. That is not the case. 

The first indication of a serious fiaw comes from the case where we do 
not put any mutations in the model. Making M^ = in (II. 2p still allows to 
perform the same analysis. In such a case, one would expect that nothing 
happens anymore at the limit as no evolution should be possible without mu- 
tations. However this is not the case and numerical simulations in particular 
show a remarquably similar behaviour between the cases with and without 
mutations (see [30] for example). This phenomenon is entirely due to the 
persistence of very small subpopulations which should actually be extinct 
but are kept by our deterministic model and can therefore re-emerge later if 
the conditions are right. 



A second important problem concerns the issue of branching. Biologically 
speaking branching is the process by which one population divides itself into 
two (and then possibly more) subpopulations. Mathematically one expects 
that at the limit, u will be a sum of moving Dirac masses 



u 



y^^ai5. 



Xi(t)- 



Branching then corresponds to the case when one Dirac mass becomes two. 
This phenomenon is one of the main motivation to study models like (II. 2p 
instead of the adaptive dynamics approach for instance (see [15] for an intro- 
duction). At the limit (II. 6p . branching occurs at infinite speed, i.e. if Dirac 
mass a5x{t) divides itself after some time to into ai5x-^(t) + Oi25x2{t) then 

— |x2(t) - Xi{t)\t=tO = +00. 

Instead probabilistic models predict a finite speed at branching... 

Those qualitative disagreements are important and should be corrected. 
One would then hope to derive models that are both able of dealing with 
very large populations and still treat correctly the small subpopulations. It 
is for the moment a completely open question of how to keep the stochastic 
effects for the small populations. But one could at least try to truncate the 
populations with less than 1 individual, which is expected to be enough to 
correct the qualitative flaws of the deterministic models. 

There are already some attempts in this direction, see [30] and very re- 
cently [29]. The proposed correction in [30] consists in studying equations 
hke 

1 1 /m7 

dtUe{t,x) = -{r{x) - I-kUe) Ue{t,x) . — + Me{Us){t, x) . (1.8) 

e £ \ Ue 

The added mortality term has the effect of killing (in times of order e) any 
population with a density less than -Ug. From the modeling point of view, 
this is quite satisfactory as it corrects most of the problems with (II. 2p . 

Unfortunately the mathematical analysis of an equation like (11.81) is for 
the moment untractable. The only situation that is understood is when Ue 
is chosen like exp{—ip/e). However this is exactly the scaling that was giving 
less than 1 individual and it is not satisfactory. Instead given our scalings, 
one would like to work with Us that are polynomial in e. 

However that would mean truncating every value less than Celoge in 
iPs, which implies at the limit, as y? < 0, everything... 



1.4 The proposed correction: Cooperative interactions 

We propose here a correction that is inspired from the case with sexual 
reproduction. In the present context of mostly asexual reproduction, it can 
however be better understood as taking into account some cooperative effects 
between the individuals. 

We can prove that it completely corrects all the abnormal behaviours of 
(11. 2p at the limit. In addition it is for the moment the only correction for 
which one can derive rigorously the limit. That means in particular that one 
can obtain numerical simulations for realistically low e. 

There exists a well known phenomenon in the case of sexual reproduction 
which drives to extinction small populations. If the population is too small 
then the probability of meeting a partner is very low and hence the birth 
rate declines. This is the idea that we follow here. 

Let us explain it first in the context of sexual reproduction. Consider a 
subpopulation with trait x, if the density of population for x, u^{x) is below 
a critical value m^ then we assume that the probability for an individual with 
trait X to meet a partner is too low. Instead this individual will reproduce 
with an individual with a different trait y such that Ue{y) is large enough. 
Typically one can expect then that y should be the closest trait to x with a 
population large enough. However as the two traits are different, the individ- 
uals are not as compatible and the corresponding birth rate should decrease 
with the distance \x — y\. 

The same phenomenon can be seen if one assumes that individuals of 
similar traits may cooperate. Selection-mutation models rightly focus on 
competition between individuals as the main interaction mechanism in order 
to observe selection (and hence evolution). Nevertheless cooperation usually 
exists as well; we assume here that it takes place at a smaller scale than the 
competitive interactions. 

We are for example lead to add to the reproduction rate a cooperative 
effect of the kind 

max f Co, /" $ (|x - y\h,) ^^ dy] , (1.9) 

for a small parameter rj^ and a decreasing $. Note that the maximum effect 
of cooperation is capped (obviously birth rates cannot blow up). 

Of course if one expects to see an effect on the small population then it is 
necessary that the reproduction rate without cooperation be negative, which 
means that the cooperation effect should be strong enough. 



It is possible to simplify fll.9p (for numerical purposes for instance) while 
keeping the same structure and the same effects. When one combines it with 
the lower value of the reproduction rate, and assume for instance a linear $, 
then it is essentially possible to replace fll.9p by a regularization of 

-Kd{x, {Ue>Us}. (1.10) 

Those are the type of corrections that we consider here, with K large enough. 

1.5 A brief overview of the various approaches for se- 
lection, mutation dynamics 

The stochastic approach is based on individual-based models. As we men- 
tioned before, they are related to evolutionary PDE models as those here 
or in [131 123] through a scaling of large population (see again |H]). Using 
a simultaneous scaling of large population and rare mutations, a stochastic 
limit process was obtained in [B] in the case of a monotype population (i.e. 
when the limit process can only be composed of a single Dirac mass), and 
in [11] when the limit population can be composed of finitely many Dirac 
masses. Other features can be added to those models, age-dependence for 
instance as in [2B] . 

At the deterministic level, one approach consists in studying a simulta- 
neous scaling of mutation and selection, in order to obtain a limit dynamics 
where transitions from a single Dirac mass to two Dirac masses could occur; 
the famous branching phenomenon that we have already mentioned. This is 
where the deterministic models presented here fit. More is said about the 
contributions in that area in the next subsection. 

Another approach consists in completely separating the two scalings. One 
then tries to directly characterize evolutionary dynamics as sums of Dirac 
masses under biologically relevant parameter scalings, instead of obtaining it 
as some limit. This is a key point in adaptive dynamics — [221 123 [2H1 [IH E]- 

A classical way of justifying this form of the solution consists in studying 
the stationary behaviour of an evolutionary model involving a scaling param- 
eter for mutations, and then letting this parameter converge to 0. The sta- 
tionary state has been proved to be composed of one or several Dirac masses 
for various models (for deterministic PDE models, see [H O [131 1231 [H] , for 
Fokker-Planck PDEs corresponding to stochastic population genetics models, 
see [3], for stochastic models, see [33], for game-theoretic models, see [12])- 



Closely related to these works are the notions of ESS (evolutionarily stable 
strategies) and CSS (convergence stable strategies) [2S1IIS], which allow one 
in some cases to characterize stable stationary states [H [131 1231 |T2] . 

In this context the phenomenon of evolutionary branching, which is es- 
pecially important for us, simply corresponds to the direct transition from 
a population composed of a single Dirac mass to a population composed of 
two Dirac masses, [28| [T9| [20] . 



1.6 An overview of the works around Eq. ( [1.21 ) 



The scalings and the first formal results have been obtained in [T7]. This 
was followed by several works on other models and on the corresponding 
Hamilton- Jacobi PDE [5l[3l]. 

The main difficulty as mentioned before is the identification of the weak 
limit of u in terms of if. There is usually one additional information which 
comes from uniform bounds on the total mass, it is that 

iaaxip{t, .) = 0, Vt. 

One solution is then to try to see the undetermined J t^t -u as a Lagrange 
multiplier for this additional constraint. 

For example if ip solves fll.6p then our constraint should imply that r—I-ku 
is non positive on the set {ip = 0} and vanishes on at least one point. 

Unfortunately, it is in general even formally not possible to identify I -ku 
with that. For instance if cp attains its maximum at one point then one has 
only constraint which is not enough. 

There are however cases of competitive interactions where this suffices. 
Suppose for instance that 1 = 1. Then I ic u is just the total mass and 
formally one simply expects that 

u{t,y) dy = max r{x). 

X s.t. ip(t,x)=0 

This is typical of so-called one resource interaction, meaning that the indi- 
viduals interact only through one average quantity. In general that means 
considering interactions of the type 

R i X, ri{x) u{t, x) dx ] , (l-H) 



with an increasing R. 

For models of the type we consider here, rigorous results (especially for 
the well posedness of the Hamilton- Jacobi eq. at the limit) mainly only exist 
in this case with just one resource, see [2] and [1], [21] (one resource but 
multidimensional traits). 

However it is possible to extend the theory, see [9] . Formally one expects 
the limit u to satisfy the following conditions 

{i) supp'u(t, .) C {v3(t, •) = 0}, 

(M)r-/*M<0, on {(/?(t, .) = 0}, (1.12) 

{iii) r — I -ku = 0, on supp u{t, .). 

This corresponds to the definition of an evolutionarily stable strategy. If one 
assumes a strong competition, i.e. that the operator I-k is positive, then 
there exists a unique measure satisfying (I1.12p (see [23] )• 

This is the approach followed here as well. With respect to the case 
with one resource, there are nevertheless additional difficulties. It is much 
harder to control the time oscillations of the reaction term, which is usually 
increasing and hence BV in time for one resource. 

2 The model studied and the results 

According to the previous considerations, we study the following equation 

dtUsit, x) = -{r-I^, Ue{t, .) - De{Ue)) Ue{t, x) + M,{Ue){t, x) , (2.1) 

where we recall that the mutation kernel M^ reads 

Me{f){x) = - [ K{z) {fix + ez) - fix)) dz, (2.2) 

for a -ft' G C^(M) such that J^zKiz) dz = 0. Following [17], one defines if^ 
as 

Ue = e"^"/^, or ipe = e logUs, (2.3) 

and obtains the equation 

dt<fe = r-I^u, + H,iy^,), (2.4) 

9 



with 

H,{f) = f K{z) (e(/(-+-^)-/(-))/- - 1) dz. (2.5) 

Jr 

First of all we need to make sure that only the traits in a compact interval are 
important (to avoid traveling waves effects for instance). This can be simply 
ensured by asking all traits out of an interval to have a negative reproduction 
rate 

3R > 0, 3ro > 0, V|x| > R, r{x) < -tq. (2.6) 

Of course the whole population should not vanish immediately, and it is 
necessary that a non negligible part be concentrated on [—R, R]. The total 
population should also be bounded which leads to 



sup / Ueit = 0, .) dx < CXD, 

' ^" ,R,, (2.7) 

3i?o iiiin ^(2^) > 0, inf / ^^(t = 0,x) dx > 0. 

We assume that the individuals interact through a strong competition 

V/ e M\R) \ {0}, / I{x - y) df{x) df{y) > 0, or / > 0. (2.8) 

This allows us to define uniquely the ESS as per 

Proposition 2.1 Assume (12. Sp . (12.61) and that r, I G C(]R). For any closed 
fi C M, there exists a unique finite nonnegative measure /i(fi) satisfying 
i) supp fi C Q 
a) r — I ic fi < in Q, r — I-k^ = Oon suppfi. 

Note however that sometimes, one may have uniqueness of the environmen- 
tal variables whereas the population measure is not unique (and (12. 8p is 
violated), most of our method would remain valid in such a case. This is 
the situation of a single resource, where almost nothing is required, see |2]. 
Nevertheless in more general situations, the conditions for which this kind of 
property holds are not currently identified. 

We need an additional assumption to make sure that the ESS can only 
be concentrated on a set of measure 

3S G C(M+) with S{0) = s.t. V/i G M\R+) 

\{x, \rix)-I^,fiix)\<u}\<Siiy). (2.9) 
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This condition is in part technical and is required to avoid some time oscil- 
lations. However it also corresponds to the natural biological idea that only 
a few traits can be present at a given time. 

It is probably hard to check (12. 9 p in specific models, but it is at least 
satisfied in large classes of parameters. One easy example is if the derivatives 
r^*^^ and I^''^ are positive (or negative) for some k. It has in fact been proved, 
see [21], that generically in r and I the ESS is discrete (a finite sum of Dirac 
masses) and hence (12. 9p should be satisfied. 

As for D^, we assume that there exists a critical scale (fe s.t. 

SUpSUp(|De| + \dxDs{t,x)\) < OO, sup sup (^e \dxxDe{t,x)\ < OO, 
£ t,x e t,x 

D,{t,x)=0 if ip,it,x)>-^, (2-10) 

\dxD,\ >2\dxr\ if d{x, {y, ip,{t,y) > -(^J) > 2(^,. 

Those assumptions are compatible with the type of corrections like (11.90 and 
(ll.lOp . They introduce a new scale in the problem —(f^ which corresponds 
to a critical population density of exp(— (^e/e). In line with our previous 
consideration, we would like to take this polynomial in e which means —(f^ ~ 
e log e or more precisely for a constant C uniform in e 

elog- <'^e<C6 log-. (2.11) 

6 e 

There is no need to be more precise on D^ to study the properties of Eq. 
(12.41) . However if one wants to identify the limit, it is necessary to specify 
what Df. should look like at the limit. Therefore we make the additional 
assumption, for any fixed / G C(M), non positive 

D, = D,(/)^min(irrf(x, {/ = 0}), Do) as e -> 0. (2.12) 

This corresponds to (ll.lOp but other shapes would be possible and would 
essentially work the same. 

Formally, we can hence expect that as e — )■ 0, Eq. (12. 4p will lead to 

9i^ = r-/^,M{(^(t,.)=0})-min(ird(a;, {/ = 0}), D^)+H{dxV). (2.13) 

where the Hamiltonian H^ became 

H{p) = f K{z) {eP' - 1) dz. (2.14) 

This is indeed what one can prove 

11 



Theorem 2.1 Assume K, r G L~ n C^{M.), J^zK{z)dz = 0, dMD on r, 
fl2.7p on the bounds for the initial population, {\2.S\i on the interaction I i\2.9\i . 
fl2.10p on Ds with (12.111) . and fl2.12p on the limit of D^, that the initial data 
Ueit = 0) > or v9£(t = 0) are C"^ , with 

mi'mi dxx^eit = ^ix) > —oo, (2-15) 

and that ips{t = 0, ■) converges to a function (f^ for the norm \\ ■ ||l°°(ir)- 

Then up to the extraction of a subsequence in e, ip/, converges to some 
continuous ({> uniformly on any compact subset of [0, T] xR and ip is a solution 
to fl2.13p almost everywhere in t, x with initial condition ip{t = 0, ■) = ^9°. In 
particular the function I -kUe converges to I -ku in L^([0,T], C(M)) for any 
p < cx), where u{t, .) = fi{{ip{t, .) = 0} is defined from ip by Prop. \2.1\ and is 
continuous in time. 

From a practical point of view, computing the solution u^ of Eq. (12. ip is 
often too costly for small e. This result allows to approximate the population 
density Ue for small e by the simpler n{{ip{t,-) = 0}), where ip may be 
obtained by a discretization of (I2.13p . in the fashion of those done in |17] . 
Rigorous numerical analysis of this kind of Hamilton- Jacobi equations is 
however still very preliminary. 

Section [3] gives a short sketch of the proof of Prop. 12.11 and can be safely 
ignored if one is familiar with [23], |32] . 

The proof of Theorem 12.11 is given in the next section and being quite 
technical is divided into several lemmas. Lemma [4. II just corresponds to the 
classical a priori estimates on equations like (12. 4p . Lemma 14.121 essentially 
follows the step of [9] and can also be skipped if the reader is familiar with 
that work. 

In the proofs below, C denotes a numerical constant which may change 
from line to line but which only depends on T, norms of the initial data or 
the coefficients, but which is always uniform on e. 

We define as usual the distance between a point and a set 

d{x, Q) = inf \x = y\. 
yen 

We also define the semi distance 6{Oi,02) between two sets Oi and O2 as 

usual by 

6{Oi,02) = sup inf \x — y\. 
xeOi y^02 
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A small S{Oi, O2) indicates that Oi is almost included in O2 and in particular 
5(Oi,02) = Oiff Oi CO2. 

We denote by M^{uj) the set of signed Radon measures on the subset u 
of M equipped with the total variation norm. 



3 Sketch of the proof of Prop. 12.1 



A complete proof of this proposition can be found in [23], and [32] . We only 
show uniqueness and give a rough skecth of the existence part. For a given 
compact set Q, assume that we can find two measures /ii and fi2 with support 
in Q and such that 

r — I -k fii = onsupp/ij, r — I -k fii < on fl. 

Compute 

> (r - I -k Hi) dfi2 = I'k{fi2- /^i) dn2. 

Summing with the symmetric term, one obtains 



/ I 'k{fi2- /ii) {dfi2 - dfii) < 0. 



However since J > then the corresponding quadratic form is positive and 
we can conclude that ^i = fi2- 

For the existence, one considers the equation 

dtfis = -{r - I k fi^) H^ + eM{us), fie{t = 0) = /i°, 

for a well chosen initial data /i°. If Q is fully discrete, Q = {xi, ...,Xn} then 
one may choose /i° = Yli ^x^- If ^ is an interval then one simply takes /i^ = 1. 
The general case is trickier and consists in taking an initial measure /i° s.t. 
for any 5 > Q and Xq € i7 

/ d^i^ > 0. 

jQnlxo-S, xo+S] 

It is then easy to obtain lower and upper bounds for the total mass. 

Note that though there is also a mutation term here, its scaling is com- 
pletely different from (12. ip . The scaling does not involve small mutations 
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and moreover the mutations are of order e (instead of 1) thus vanishing at 
the hmit. 

Passing to the weak hmit, /i^ — >■ /i, gives a measure that has the desired 
properties (but is by no means easy to show). 



4 Proof of Theorem 12.1 



4.1 A Priori estimates 

We start by stating and proving the obvious a priori estimates that one can 
obtain for the problem. Those essentially follow the lines of previous works. 
We show the following estimates on the solution to (12. 4p 



Lemma 4.1 Let </?£ he a solution to (12. 4p with the assumptions of Theorem 
\2.1\ Then for any T > 

||<9tV^£||L°°([0,T]xR) + ||<9a;V5£||L°°{[0,T],L°°(R)) < Ct, (4.1) 

Vt<T, xgM, 9,,^,(t,x)>-£, H,{^,)>-Cf^, (4.2) 

Vt <T, ^ < /jgM^(t,a;)rfx < C, i^^{t, x) < e\ogl/e + C e, (4.3) 

where C only depends on the time T, J^e"^^^^^^^ dx, ||c?xV^el|L°°(R) ond the 
infimum of d^xfei^) ■ Finally ip^ has level sets, uniformly bounded in e. 

Proof. Almost all proofs here are taken directly from [9] and reproduced 
for the sake of completeness. Some resulting bounds are different and much 
worse than in this former article, more precisely the lower bounds on dxx^e 
and H^. Nevertheless even for those, the proofs, i.e. the way to obtain the 
bounds, are very close. As such we may skip some technical details. 

Step 0: Upper Bound on the total mass. First notice that because of 
([22]), there exists i? > s.t. 

V|a;| > i?, r{x) — I -k u^it, x) < —tq. 

Let tp he a, smooth test function with support in |x| > R, taking values in 
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[0, 1] and equal to 1 on |x| > i? + 1. Using the previous bound, we compute 

— / il){x)u^{t^x) dx < / ijj{x)ue{t,x) dx 

dt JiR s J^ 

1 [ 

H — / K{z){il){x — ez) ~ il){x))ue{t,x) dzdx 

^ Jm? 

< / ip{x)ui;{t,x) dx + C / Uf.{t,x)dx. 

^ Jr Jr 

On the other hand, on the bounded domain [—R — 1, -R + 1] as I{x) > 
for all X, one has for some constant C 

V|x|<i?+1, Ii.Ue{t,x)>C I {l-^)Uedx. 

JR 

Therefore with the same kind of estimate 
^/(l-*W)«.(,x)d.<c/,..(,x)dx 

dt Its Jm 



H — sup r—C {1 — ^) Ue dx \ ( 
^ \ Jr J Jr 



[1 — ip{x)) Us{t, x) dx. 
fc \ ^R / Jr 

Summing the two 

dt Jr £ \ Jr J Jr 

/ ipix) Us{t, x) dx + C / Us{t,x)dx. 

^ Jr Jr 

Since the sum of the first two terms of the r.h.s. is negative if J Ue is larger 
than a constant independent of e, this shows that J Us{t,x) dx remains uni- 
formly bounded on any finite time interval. 

Step 1: Bound on dx^e- This is a classical bound for solutions to Hamilton- 
Jacobi equations and we have to check that it remains true uniformly at the 
e level. We follow exactly [9] for instance. Compute 

dtdxfe = dxT — dxl ^ «£ — dxDg 

^,[t,a^+ez)-^,(t,^) dx(Pe(t, x + ez) - dxiPeJt, x) (4.4) 
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By our assumptions and the upper bound on the total mass 

\dxr — dxl * Me — dxDe\ < C. 

First note that this shows that ||(9j^(y9e(t, .)IU°° remains finite over a (possibly 
very short) time interval [0, tg]. 

Now we use the classical maximum principle. Fix t G [0,T] such that 
Ce,t ■= \\dx'^e{t,-)\\L°°iu) < oo. For any x e M such that dx(Pe(t,x) > 
supydxfeitjy) — a, where the constant a > will be specified later, we 
have 

dtdxipeit,x) <C+ [ K{z)e\'\^''^-dz<c(l + -e^^'A. 
Therefore, choosing a = ee"*" *"*'=, we obtain 

— sup dx(Ps{t,x) < C, 
dt X 

for a constant C independent of t < tg and of e. Using a similar argument for 
the minimum, we deduce that t^ > T and that dx'-Pe is bounded on [0, T] x M 
by a constant depending only on T and ||9a;V9°||i^oo(iK). 

Step 2: First hound on Hg{ips) and bounds on dt^Pe and (pe- First remark 
that 

- f K{z) dz < H,{ip,{t)){x) = I K{z) e '^-'*--+",'-^-<'--> rfz - /" K{z) dz 

Jr Jr Jr 

< / K{z) el^l I19-¥'.IIl-([o,t],e) clz < c. 



This is not optimal, more precisely the lower bound is atrocious, but will 
suffice for the moment. 
Directly from Eq. fl2^ . 

\9t^s\ < sup \r\ + sup |/| / Usdx + \De\ + C < C, 

Jr 

hence ending the proof of the full Lipschitz bound on cp^. 

To get the upper bound on ip^, we use this Lipschitz bound to get 

^e{t,y) > (Peit,x) -C\y-x\, 
16 



so 

f Ueit, y)dy> I e^^(*'^)/" e"^ l^"^!/" dy > 2C-^ e e^-(*'^)/^ 

Hence the bound on the total mass yields that ^e < e log 1/e + C e. 

Step 3: Lower hounds on d^x^e C'l^d H^. Here we start paying for the 
introduction of an additional mortality term as the second derivative of D^ 
is not bounded uniformly in e. As before we use a maximum principle, from 
([22D and (1210|) 

C f ^.. . ^.(t,^+ez)-^,(t.x) d^^ipAt.x + ez) - d^^ipAt.x) 
OtOxx^e> -^ + I K{z)e = dz 



(f,^+Ez)-ye(t,:r) ((9a;(y9j(t, X + ez) — d^feit, x)) 



2 



+ / K{z)e ^ dz. 

Jr £ 

The last term is of course non negative and so with the same argument as 
before, we get 

— midxx^eit,x) >- — . 
at X ip^ 

This proves the uniform lower bound on d^x^e- Let us turn to the sharp 
lower bound on H^{(f^). Let us write 

Hsi^e) > / K{z) exp ( / zdx^e{t,x + eze)d6] dz - K{z)dz. 
Jr \Jo J Jr 

By differentiating once more 

/ z dx(ps{t,x + 6ze) d6 > zdx(peit,x) 
Jo 

+ z ezedxxV>e{t,x + e'eze)de'de 
Jo Jo 

> zdxipe(t,x) -C ^z'^. 
Eventually 

Hei^s) > / K{z)exp{zdx(Pe{t,x) -Cez'^/<fs)dz- / K{z)dz 



>H{dx<fe{t,x))-C—, 
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where H is defined as in fl2.14p and since K is compactly supported. Because 
we assumed that j^ zK{z) dz = 0, we have H{p) > for any p, which gives 
the final bound. 

Step 4- ^e has uniform compact level sets. 
Observe that ipsit = 0, x) — )■ — oo when x — )■ ±oo since J^u^it = 0, x) dx < 
oo and dx^Peit = 0) is bounded. Because of the uniform convergence of 
ipe{t = 0) to ip^ on M, one moreover deduces that the convergence ips{t = 
0,x) — 7- — oo is uniform in e. 

Since dx^Pe £ -^°°([0,7'],K), this remains true at any time t G [0, T], 
uniformly in e. 

Therefore, the set 

Qe ■■= {{t, x) G [0, T] X M : ip,{t, x) > -1} 

is bounded for any /, uniformly in e. 

Step 5: Lower bound on the total mass and ips- 
By the previous steps we know that 

M,{u,)>^H,>-C^. 

Note that 0^ = whenever u^ > e^'^^'^. Therefore by (12.71) . choose Rq 
s.t. mm\x\<Kr > and integrate (12. ip 

dt J\x\<K ^ \\^\<K 'Pe Jk J J\x\<K 

This implies that by (12. 7p 
log / u,St, .) <,. > -C + i /' f mm . - C ^ - C / ..(.,.) i.) is, 

J\x\<K ^ Jo \\^\<K "Pe Jr J 

or obviously 

log / u,(t, .) ,. > _C + i /'( mm . - C ^ - C / ..(». ^) iA *. 

Jr £ Jo \\^\<K ipe Jr J 

This allows us to conclude that the total mass remains bounded from below 
uniformly in e. 

In particular by step 4, this lower bound means that max^^^ > —ips- 
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4.2 Passing to the limit in the equation: First steps 



By the uniform bounds provided by Lemma HTTj we can extract a subsequence 
in e (still denoted with e), and find a function yj on [0,T] x M such that 
9t(^ e L°°([0,T] X M), 9^^ G L°°([0,T] x M), with max(^ = satisfying by 
the Arzela-Ascoli theorem 



^e > ^ 



uniformly in C{K) for any compact K of [0, T] x M, (4.5) 



Note that by the upper and lower bounds on yj^, one has maxip = 0. 

Since u^ is uniformly bounded in L^, it converges (still after an extraction) 
in the weak-* topology of measures to some u G -Z^°°([0, T], M^(M)). This 
in turn implies a weak convergence of I -k u^. Note that it is indeed only a 
weak convergence in spite of the convolution because it regularizes only in x 
and time oscillations are still possible. 

Using the notion of viscosity solution (see [2] for instance) and the uniform 
bound on d^ip;;, we may obtain the convergence of H^{ip^) to H^dr^ip). 

Finally as D^ is uniformly bounded and can hence be assumed to converge 
weakly to some D G L°°, we obtain 

dtip = r-I^ u{t, .)-D + H{d,ip). (4.6) 

Unfortunately D and u are still unidentified. This is now our aim but it will 
require a much more precise understanding and control of the time oscilla- 
tions of the set where the population is concentrated. 



4.3 Continuity in time of the set {x | (/? = 0} 

At the limit, it is possible to show that the points where the population is 
concentrated move at most at a finite speed given by 

l^ = 2 sup sup -— e''^K{z)dz-l. 

£ l€l<2||<9:,</3^||iOO Kl J 

We comment on that again at the very end of the proof. However we do not 
have yet enough tools to prove such a strong statement. For the time being 
we will be satisfied in proving a continuity result. The first step is to do it 
at the limit 
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Lemma 4.2 Assume that for some point Xq, som,e t > to and S > 

d{xo, {ip{s,.) = 0})>6, Vse [to, t], 

then {pit, Xq) < 0. 

Before turning to the next result, let us point out that Lemma 1121 actually 
implies the following 

3r G C{R+) with r(0) = 0, si. Vs > t, 

Vx G {(/?(s, .) = 0}, 3y e{<^{t,.) =0} with \y-x\<T{s -t). ^'' 

One of the main idea in the following proofs is to follow the character- 
istics corresponding to the Hamiltonian H. However Eq. (12 ■4p involves the 
modified Hamiltonian if^ which does not have characteristics per se (it is 
non local for example). Some modifications are hence needed and follow the 
usual ideas for parabolic problems, which is why we use the following lemma 

Lemma 4.3 Introduce the intermediary scale a/F|^7| = ^/^ = ^'^e ^^f^d 
consider any interval [a, h]. Define a{t) = a + {t — t^) V/2, b{t) = b — {t — 
to) V/2. We denote 

mJt) = max u:>Jt,^. 

Then 

ms{t) - m^{to) > / ( min {r{x) - I -k^Ue{s,x) - D^) (fie) ds 



nieit) < max v^£(to, .) 



/tg V [a(t), h{t)] 

+ / max (r(x) — I -k^Ug^t^x) — Di; + 2e) ds. 

Jto Hs)-£V'e, b{s)+eips] 

Proof of Lemma 14.31 The part 

min {r{x) - I -k^ Ue{s, x) - D^) 1^^ < —me{s) 

[a{t), fe(t)] s at 

is direct once one recalls that for any x G [a{t), b{t)] then 

dt(fe{s,x) = r{x) - I -k^Ue{s,x) - De + Hs{fe){s,x) 

C 

> min (r(x) — I -k^uJs^x) — D^) (^e, 
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(4.9) 



by the lower bound fl4.2p . 

Let us turn now to the upper bound which is trickier and involves the 
velocity V. 

Consider Xe a regularization of |||9kV^£||l°° () + , with ()+ the positive part, 
at the intermediary scale eips , i-e. 

Xe{x) = 0, if X < 0, 

3 __ (4.8) 

together with the general bounds 

3 3 

C 

< ^^X) < 2 \\d:rVs\\L^^x>0 + e<^s, \Xe{^)\ < — ^• 

With that one defines 

Ipeit, X) = -Xeiait) -X)- Xeix - b{t)). 

Note that ipe satisfies the following inequation 

V 

This has for consequence that Vz G supp-ft' C [— 1, 1] 
/ exp (mi^^ + e^)-Ui^^)\\ K^,^ ^, _ 1 

< / e"^^'^^^*'^') K{z) dz~l + e<— |(9^Ve| + e. 

As ipe and iIj^ are smooth functions, it is easy to study 

fheit) = snp{ip,{t, .) + ^^(t, .)) = max{(p^{t, .) + ^^(t, .)). 

In general this maximum is attained at one (or several) point x G uJm{t)- 
First note that at such a point, one has 

d^{ips{t,x) +ilJs{t,x)) = 0. 
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By fl4.8p this implies that 

cu^t) C [a{t) - e^„ b{t) + e^e]. (4.11) 

Now in general 

—m,{t)< sup dt{ipsit,x) + ilJe{t,x)). 
Note that by the definition of the maximum, Vz G supp K 

iPe{t, X + ez) < ipe{t, X) + Ipeit, x) - Ipeit, X + E z) . 

Now for X G ujm{r) 

HA^.){t, X) = I exp ( ■^'"■" + "^>-'^'"-"> ) A-(.) d. - 1 
< / exp (-Mt.n: + ez) + Mt,x)\ ^^^^^ ^^ _ ^ 



V _ 



Therefore 

d 



me{t) < sup {r{x)- I -k^Ue{t,x)- D^)+e 

"* X&UJ„,{t) 

< sup (r(x) -/7kra,M£(t,x) -De) +£, 

xG[a{t)-Ce(^£, b{t)+Ce(pe] 

< sup (r(x) — 1*2, Me(t, x) — De) + 2e 

a::G[a{t), b{t)\ 



bygHD. 

Note that by ^Ei and fimjl 



me(t) < me{t) < max __ me{t), 



which allows to conclude. □ 

Let us now turn to Lemma 14.21 
Proof of Lemma 14.21 
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Step 1: Semi- continuity of {if = 0} 
Let us start with the following crucial observation 

Vt, 3Tt e C(M+) with ri(0) = 0, s.t. Vs > t, 

Vx G {^(s, .) = 0}, 3?/ G {y?(t, .) = 0} with ly - x| < n(s - t). 

This is a sort of semi-continuity for {(p = 0}. It is proved very simply by 
contradiction. If it were not true, then 

3t, 3ro > 0, 3s„ -^ t, s„ > t, 3a;„ G {ip{sn, •) = 0}, 
rf(x„, {v9(t, .) =0}) > To, 

where d{x,co) = infygi^ \x — y\ is the usual distance. 

Since all the x„ belong to a compact set , we can extract a converging 
subsequence Xn -^ x. As </? is continuous, v^(t,x) = or a; G {(/^(t, •) = 0}. 
On the other hand one would also have d{x, {v?(t, .) = 0}) > tq which is 
contradictory. 

Therefore the result of the lemma is obviously true if t is such that Tt^, (t) < 
6. 

Step 2: The connection between Tt and the Lemma. 
One does not have in general a uniform control on tiq. In that case, one 
would find a sequence t„, a number 6 > s.t. V?7 > 0, lim sup r^^^ (77) > S. 

The result of Lemma 14.21 would precisely rule this out and we argue by 
contradiction. Denote by t the first time when such a jump occurs. That 
means that 

• T can be chosen uniform till t, i.e. 3r s.t. Vs < t and V77 G [0, t — s], 

• There is a jump at t of size 6 > 0, which means 3xo, 3to < t s.t. 
ip{t,xo) = but 

d{xo, {ip{s,.) = 0})>6, VsG [to, t]. 

Note that one can take to as close to t as one wishes and in particular we 
may freely assume that t — to is small enough s.t. T{t — to) < S/8. 

Step 3: The contradiction. 
First of all, note that as y^e — > V^ in L°° norm, one has that for e small enough 
and some interval Jg with iJgl — )■ as £ — ;■ 

d{xQ, {v?e(s, .) > -ipe}) > 2' ^■'^ ^ [^0, t] \ le. 
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Indeed one observes that 
And by the continuity of (p 

HWi-'^^ ■) > -Ve - Wve - v^IIl-}, {v5(s, .) = 0}) — > 0, 

for almost every s. 

In particular we point out that one may define le with \I^\ — )■ and 
r^g -> s.t. 

5{{^e{s, .) > -^e}, Wis, .) = 0}) < r],, Vs e [to, t] \ I,. (4.13) 

Now define the following interval 

his) = [xo - 6 /A + is- to) V/2, xo + 6/A + (s - to) V/2]. 

Moreover denote by yo the closest point on the left from xq in {y9(to, .) = 0}. 
Assume for instance that yo < xq. 

Next for any e small enough, and any s G [to, t], denote 

yeis) = sup{a; < xq, ^eis,x) > -(^J. 

By (14.1 31) . one has that ?/e(s) < yo + r]^ + r(s — to) for any s G [to, t] \ Jg. 
The same result can be obtained for the closest point on the right. 

By the assumption on xq in the Lemma, this implies that Ji remains at 
a distance larger than 6/4 > C ipe oi {peis, .) > —p>e}- 

By the properties of De, we deduce that for any s G [so, t], 

maxr-JTkr^. ^^(s,.) - D^ < r(|/e) - I -k^ Ueis,ye) - Deis, ye) -C 5/A 

h{s) 

< riyo) - I-k^ Ueis, yo) - Deis, yo) - C (5/4 - rje - r(s - to)) 

< riyo) - I -k^ Ueis, yo) - Deis, yo) - C 5/16, 

by choosing e small enough. 

Now let us apply Lemma H75] to Ji to deduce that 

maxv3£(t, .)< / iriyo) - I'k^Ueis,yo) - Deis,yo))ds 

^W J to 

-C — it-to) + C\Ie\+ max ipeito,-)- 

16 [xo-35/4, X0+3S/4] 
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Note that by K^ 

Ce 
dtife{s,yo) >r{yo) - I 'kxUs{s,yo) - Ds{s,yo) - -^. 

So 

^it,Xo) < ||V9- V3e|U°o +maXV3e(t, .) <'fe{t,yo) -<fe{to,yo) ^ C — {t - to) 

hit) lb 

+C\Ie\ + \\ip-iPe\\L°-+ max '-Peito,.). 

[X0-3S/4, a:o+3<5/4] 

To conclude, just observe that ip{tQ,yo) = so that 

'feit,yo)~^e{to,yo)+C\Is\ < \\^-^e\\L-°+£ \'^oge\+C \Is\ — ^ 0, as s ^ 0. 
Therefore taking e small enough, one concludes that 

V?(t, xo) < max ip(to, .) + \\(p - (pAlL'^ < 0, 

N-35/4, xo+35/4] 

which gives the desired contradiction. □ 

4.4 Control of the oscillations 

Before obtaining the continuity of the set {ips > —(fs} which is our main 
goal, we have to control the oscillations in time of the reaction term. 

Lemma 4.4 There exists a positive function ae(s) with J^ ^^(s) ds < 1, s.t. 
for any to ^ [0, T] and any t > to 

1 /•* 

\\I -kUeit,.) - I -kUe(to,.)\\L^ < —= / ^^(s) rfs. 

V ^ Jto 

With this control, we can first show that the reaction term r—I-ku^ essentially 
vanishes on the set {v?e > —^e] 

Lemma 4.5 There exists e ^^ 0, Is composed of at most C e~^^^ intervals of 
size e^l^ s.t. for any to G [0, T\\Is, for any xo s.t. 

3t G [to, to + e^l\ <Pe{t, Xo) > -<Pe 

then one has for any s G [to, to + s^^"^], 

r(xo) - I-kUe{s,Xo) < e. 
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This last lemma enables us to be very precise concerning the growth rate at 
any point without mutations 

Lemma 4.6 For any to and t > to with t — to > C e, for any xo and any 
u > C Vfc then 

either / (r(xo) — / * (s, Xq) — -De(s, Xq)) ds < —C{t — to) {v"^ — e), 

Jto 



/ De{s,xo)ds<C{t~to)v. 

Jtn 



or 

'to 



With this control on the growth rate, it is possible to make sure that at most 
times the set {^p{t, .) = 0} is of measure or more precisely 

Lemma 4.7 For any to G [0, T], define the set Vtto by xo G Qto iff 

Xo e {^{to, .) = 0}, 3t„ -^ t, XnG {(fitn, •) = 0} S.t. X^ -> Xq. 

Then \Qt\ = 0. 

We now turn to the proofs of those results 
Proof of Lemma 14.41 The proof uses one of the entropies of the system 
(and a different one from [9j). Start by the following 



d_ 
di 



Ue{t, x) I{x — y) Ue{t, y) dx dy — 2 r Uf.{t, x) dx 
- {r ~ I -k Ueit, .)) Ueit, x) (/ -k {ueit, .) - r) dx 



2 f 

H — / D^(t, x) Ueit, x) (/ -k Ueit, .) — r) dx 

+ I MeiUe)it,x)iI^Ueit,.)-r)dx. 



Note that as D^it, x) = if yje > —ipe then 

DeUe<Ce-^^'' <Ce. 
Similarly for any smooth function ?/;, by a simple change of variable 

MeiUe)^ix) dx = - I Kiz)Ueit,x) (V^(x - £ z) - V^(x)) dx < C £ || V^ ||loo . 
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Hence one first obtains that 



£ JO 
Now note that 



1 r^ , 

{r-Ii. Ue{t, .)f Ue{t, x) dx dt < C. (AAA) 



1 /•* 

/ -k Ue{t, .) - I-k Ue(to, ■) =- I -k{{r - I -k Ue{s, .)) Ue{s, .)) ds 

^ J to 

1 /•* 

- - / I-k{De{s,.)Ue{s,.))ds 

^ Jto 

+ / I^Me{Ue{s,.))ds. 

With the same bounds as before, one obtains that 

\\I -k{Ueit) -Ue{to))\\Loo < C (t-to) + - / / \r - I -k Ue{s, .)\ U^i-S, .) dx ds 

^ Jto J 

<Cit-to) + ^ I Q fir-I^Ue{s,.)fUe{s,.)dx 

by Caucliy-Schwartz. Defining 

1/2 



ae = C + i- {r - I-k Us{s, .))^ Us{s, .) dx j 



) 



allows to conclude. □ 

Proof of Lemma 14.51 First of all, consider the intervals Jj = [2i e^/^, 2(i + 

1) e^/^] for i = . . .T e^^^^ /2. As a^ is uniformly bounded in L^ then there 
are at most C e~^^^ indices i s.t. 

xl{s) ds > £-3/^ ( f «,(S) ds^ > £3/^ 

that is 

^ f aeis) ds > e^/^^ 
V ^ Ju 

Define h as the union of the intervals h = \l%e^l^-2e^l^, 2(i + l) e^l^^2e^l^\ 
for such indices i. 
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Now for any tg ^ h ^'Hd any t E [to, to + £^^^] , one has that 

e Jt-£3/i 

Therefore by Lemma S31 I-kUei-s) changes by at most e^^^^ over the interval 

[t-e^/\ t + e'/% 
Define 

Choosing any xq s.t. v^e(t, Xq) > — <^e and any s G [t — e^^^, t + e^^''], we 
argue by contradiction. If 

r(xo) - I-kUe{s,XQ) > e. 

Then at any other s' G [t — e^''^, t + e'^^^], and in particular for s' = to, 

r(xo) - / * M£(s', xo) > C — + e^/^ 
Note that D^lt, xq) = and hence 

dtipsit, Xq) > r{xo) - I-k Ue{to, Xq) - C — > E^^^. 

By continuity that means that ips{s',xo) is increasing in s' over a (possibly 
very small) interval [t, t + u]. Define z/ such as to make this interval maximal. 
At t + u, ip^it + z/, Xq) > v^e(t, Xq) > —^e which means that D^{t + z/, xq) = 
and if t + z/ < t + e^/"^ 

5tV5£(t + z/,xo) > r(xo) -/*M£(t + z/,xo) -C— > e^^^. 

This implies that t + u > t + e^/^, and hence that (/9e(s',Xo) > — v^£ on any 
s' e [t, t + £:3/4]. Therefore 

9i(/5.(s',xo) >e^'\ 
or finally 

¥P,(t + £=^^Xo)>-(^. + £'/', 

which is impossible as the maximum of y^^ is of order e \ log e\. □ 
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We turn to the next proof 
Proof of Lemma 14.61 For any s G [to, t] \ I^, one has either that Xq is at 
distance less than u of {fe{s, ■) > —<^e} or it is at distance larger than i/. 
Let us decompose accordingly [to, t] into Jg fl Ji, fl Ji, where I^ consists of the 
set where xq is at distance less than u of {(Pe{s, .) > — <^e}. 

We observe that on J^, by the properties of D^ and since on {(Pe^s, .) > 
— (^e}, one has by Lemma [4.51 that r — / -ku^is, .) is less than e then 

r(Xo) — I -k Me(s, Xq) — D^{s, Xq) < S — C U < —C V. 

On the other hand, on ly then still by 14. 5[ one has that 

r(xo) -Ik Ue{s, Xo) - Ds{s, Xq) < E, 

and moreover 

De{s,xo) < Cu. 

Therefore 

{r{xo) — I 'kUe{s,Xo) — Di;{s, Xo)) ds < \Iir\ + C e \Iy\ — C u \Ji,\, 



'to 

and 

/ D,{s,Xo)<Ciym+C\J,\. 

Jto 

We recall that \Iir\ < e^^^ which is asymptotically much smaller than e^. 
Therefore either \Jy\ > {t — to) u and we find the first possibility or \Jy\ < 
(t — to) z/ and we find the second one. □ 

The next proof is 
Proof of Lemma 14.71 Note first of all that by Lemma Hl2l one may in fact 
take any sequence t„ — )■ to for any xq G flto- 

Choose any to, any r] > and any xq in {ip(to, .) = 0}. Observe that for 
any t > to 

^e{t,Xo) >(Pe{to,Xo) + / {r - I k Ue{s, Xo) - Ds{s, Xq)) ds - C ^ . 

Jto ^e 

Therefore we first deduce that 

* e 

{r - 1 -kUeis, Xo) - De{s, Xo)) ds < C e\ loge] +C (t-to) — + \\^e-V\\Loo- 
to "f^ 

(4.15) 
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Next find Xt the closest point from xq s.t. <f{xt) = 0. By definition 
Zj — Xo — i- 0. We define the small interval around Xt for tg < s <t 

Is = [xt - (t - s)V/2, xt + {t- s)V/2]. 

By Lemma 14.31 

-\\v~Ve\\L°°<^\'^oge\+ / {r{xo) - I-k^Us{s,Xo))ds 

Jto 

+ C{e+\xo-xt\ + {t~ to)) (t - to). 
This implies that 

/ (r(xo) - I-k^ Us{s, Xq)) ds>-\\ip-ipe\\L^-e\ loge| 

Jto 

-C{e+\xo- xt\ + (t - to)) (t - to). 

(4.16) 

Now define z/ s.t. 

Ciu^-e)=^^^^^ + '-j^ + Cie-+\xo-x,\ + it-to)). 
t — to t — to 

Note that as long as e is small enough one has that u > e. Moreover u can be 
chosen arbitrarily small by taking t close enough to to. By applying Lemma 
14.61 for this z/, fl4.16p means that we are necessarily in the second case which 
means 

f D,{s,Xo)ds<C{t-to)u. 

Jto 

We combine this inequality to fl4.15p to deduce that 

I {r-I^ Ue{s, xo)) ds<C{t- to) (— + z/' j . (4.17) 

Now, we define the average 

1 /■* 

^le = — / Us{s,Xo)ds 

and from (I4.16P and (I4.17p . we get 

|r(xo)-/*/i£(xo)| <C ( — + Z/M , Vxo G^o- 
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By the fundamental property (12.91) . we deduce that 

\ntJ<S{C (— + z/' 



To conclude take t ^ Iq s.t. z/ — )■ and then e small enough to obtain that 
|fitol = 0. n 

4.5 Continuity in time of the set {x [(yi^g > —(pe) 

As a first consequence of our control over the oscillations of the reaction term, 
one may simply obtain the continuity in time of the set {ipe > ~'^e} 

Lemma 4.8 There exists Is with |/e| — )■ as e — )■ s.t. if for some to € 
[0, T] \ Jg, any point Xq, any t > t^ and 5 > 

d(xo, {Ve{s, ■) > -<^J) >5, V s G [to, t], 
then for e small enough with respect to 6, 

iPe{t, Xq) < -(ps -C6{t -to). 

for some constant C independent of e. 

Unfortunately Lemma 14.81 does not imply a result like (14. 7p with a function 
T uniform in e. It guarantees that there cannot be jumps of significant size 
in the set {ips > —'fe} but that set could still be propagated very fast. 

However we can combine it with Lemma 14.71 to finally deduce a uniform 
(in e) continuity in time for the support of {(pe > ~'fe}- 

Lemma 4.9 There exists I^ with \Is\ — )■ as e — )■ s.t. for any to G 
[0, T] \ Je, W6 > 0, there exists r > s.t. Vxq, if 

d{xo, {'Peito,-) > -'^e}) > 5, 
then for e small enough with respect to 6 and r, for any t < to + t 

iPeit, Xo) < -^e - C (t - to) < -ipe- 



Lemma 14.91 in particular directly implies 

31 



Corollary 4.10 For any to E [0, T] \ I,, Bn, E C([0, T]) with ^^(O) = 
s.t. for any t > to 



We turn to the proofs 
Proof of Lemma 14.81 Let first characterize the "good" to- We define 

Je = {to, 3t>0with |/,n[to, t]\ >£|t-to|}. (4.18) 

Note that of course I^ C Is (just take t = to + e^^^). We can estimate |/e| the 
following way: If to E le then there exists an index fc > s.t. to E I^ where 
I^ is defined as the union of the intervals [i 2'' e^/^, (? + 1) 2*^ e^^^] for those i 
s.t. 

I J, n [z 2' s'"\ % 2"^ e""^ + 2^ £3/4] I > 2'=-3 ^^3/4_ 

Each /^ is composed of at most rik such intervals and since \Is\ < C e'^^^^ one 
has that 

nu<2^-''Ce-^"'e-\ 

Consequently 

\i^,\<Ce^"'e-\ 

Since k < C \ logej, we conclude that 

\ie\ <Ce^'^e-^\\ogel 

which indeed converges to as e — )■ 0. 

Now define y^{s) similarly as in the proof of Lemma [4.21 by 

y^{s) = sup{y < xo, ^e{s,y) > -ipe}, 

and y^ by 

y^{s) = sup{y > xo, ^e{s,y) > -(^J. 

Define again the interval 

his) = [xo - 5/2 + {s- to) V/2, xo + 5/2 + {s - to) V/2]. 
It remains non empty for s E [to, t]. 
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For s G [to, t]\ hi one has that 

mii^{r -I-k^Ue{s,.) -D^) < maxr - I -k^u^^s, .) - Ds{s,.)\ + 
ii{s) ± 'y= ('^' 

-CS. 
Apply now Lemma I4.5[ if s ^ /g then this imphes that 

max(r — / -k^ We(s, .) — D^) < e — C 6 < —C6, 

his) 

provided that e is small enough with respect to 6. 
We now apply Lemma 1^3) to Ji to deduce that 

max(^,(t,.) <2eit - to) - C 6 (t - to - \I,n[to, t]\) 

hit) 

+ C\I,n[to, t]\+ max if, (to,.). 

[X0-3S/4:, xo+35/4] 

Now since to ^ h then 

\ie n [to, t]\ < sit - to), t - to - 14 n [to, t]\ > ^-^. 
Hence 

maxifeit, .) < -^e - (t - to){C5 - C e) < -^e - C 5 {t - to), 

hit) 

provided that e is small enough with respect to 6, which concludes the proof. 

D 

Proof of Lemma 14.91 By Lemma 14.21 and Lemma 14. 7[ we know that for 
a.e. to then 

\Wito, ■)}\ = 0. 

Let us choose such a to with as well to ^ I, and denote the corresponding 
interval I,. Now for any 6, divide the domain in intervals Jj of size 6/2. And 
denote 

m = supinf (y9(to, .). 

i ^i 

By the previous property m < 0. Define r = — m/(3||9?||vi/i.°°)- 

Now for any xq at distance larger than 5 from {y9e(to, .) > — f^e}, there 
exists Hi < xo and 1/2 > xo s.t. 

\yi-xo\>5/2, f{to,yi)<m. 
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Obviously for any t < Iq + r, 

(p{t,yi) < 2m/3, 
and hence for any y s.t. \y — yi\ < t then 

ip{t,y) <m/3, (Pe{t,y) <m/3 + \\(p-ipe\\L°° < -'^e, (4.19) 

provided e is chosen small enough. 

Now we apply a first time Lemma 14.81 with 6 = t which means that 
for e small enough with respect to r, there cannot be a jump of size r in 
{(fe > ~'fe}- If for some t <tQ + t 

bi, 2/2] n {if.it, .) > -(^j 7^ 0, 

then by Lemma [4. 81 with 6 = r, it means that for some s G [to, t], one can 
find a point y with \y — yi\ < r or ||/ — 1/2I < t s.t. 

which contradicts fl4.19p . 

Hence we know that xq stays at distance S/2 of the set {ips > —Ve}- We 
apply once more Lemma 14.81 to conclude. □ 

4.6 Passing to the limit in the equation: Final steps 

The first step is to characterize the limit of D;,. Of course there is a natural 
candidate for that which is 

D = mm{Kd{x, {<^ = 0}), Dq). (4.20) 

By fl4.13p . we essentially know that at most times t, the set {y^eit, .) > — <^e} 
is included in a neighborhood of {ip(t, .) = 0}. However to prove fl4.20p . 
we need to show the contrary namely that {(p{t, .) = 0} is included in a 
neighborhood of {feit, .) > —ife} at most times t. This is what the following 
lemma shows 

Lemma 4.11 For any fixed S > 0, denote 

Je = {te [0, T], 6{{v{t, .) = 0}, {Ve{t, .) > -Ve}) > 6}. 

Then J^ — ?■ as e -^ and consequently fl4.20p holds. 
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The last step is to identify the weak hmit u oi u^. Fortunately we now have 
all the required tools and we can follow the same ideas as in f9j. Therefore 
let us define 

/ii = ^(t)=/i({<^(t,.) = 0}), 

where /x is given by Prop. 12.11 

We prove the intermediary result 

Lemma 4.12 For any fixed t^ G [0, T]\/£, there exists a function a G C(M+) 
with (t(0) = s.t. 

/ / \I •ku£{s,x) — I -k fi(tQ,x)\'^ dx ds < C {t — to)T{t — to) + a{£), 
J to JR 

where t is given by (14.71) . 

Then simply by passing to the limit one deduces that 

Corollary 4.13 For any fixed to, 

/ / \I •ku{s,x) — I -k fi(to,x)\'^ dxds < {t — to)T{t — to). 
J to Jm. 



This means that any to is a Lebesgue point oi I-ku and then necessarily that 

I^u(to,.) = I'kfi{to,.), 

which as / > lets us conclude that u = fi. We have identified all the terms 
in the limiting equation which reads 

dtip = r-I^fi{{ip{t,.) = 0})-mm{Kd{x, {^ = 0}), Do)+H{d,ip), (4.21) 

which finishes the proof of our main theorem. It only remains to give the 
proofs of Lemmas 14.111 and 14.121 

Proof of Lemma 14.111 Choose a time to ^ [0, T] \ Is- Apply Corollary 
14. 101 and Lemma [4.91 and find t s.t. r^,, (t — to) < S/2 and the r from Lemma 
14.91 satisfies t < to + r. 

Then for any s G [to, t] and any Xq s.t. 

dixo, Weis,.) > -^e}) > 5. 
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First by Corollary KM 

d{xo, {fsito,-) > -<^J) >6/2. 
One has then by Lemma [4.91 that on [to, t], 

iPe{s,Xo) < -ipe -C6{S -to). 

That implies that on s G [to + C 6^^ \\(p — (Pi;\\l°°, t] then 

xo^ {v{s,.) = 0}. 
Therefore for any s E [to + C 6^^ \\lp — (P£\\l°°, t] 

6{{^{S,.) = 0}, {^e{s,.)>-^e})<S. 

We conclude that for any to ^ le there exists t > chosen uniformly in e s.t. 

Jen[to, t] C [to, to + C6'^\\ip-ipe\\L^]. 

which implies that | J^l — )■ 0. □ 

Proof of Lemma 14.121 The steps are mostly the same as in |!9]. We give 
them for the sake of completeness but we do not repeat all the details here. 
Part of the proof can also be simplified by using the more precise results that 
were proved before in the present article. 

Step 1: The functional. 
We look at the evolution of 

F^{t) = / logUe{t,x)dHto{x) = - lfe{t,x)dfit,Xx), 

Jr ^ JR 

for s > to. Compute 

—Feit) =- / {r- I-kUe~ D,) dfitoix) 

dt e J^ 

+ - / Hs{(fe{t))UsdHto- 

^ jR 



Write 



- / {r - I -kUe - De) dutoix) = -n Ue{t,x)dx 

£ Jr dt J^ 



{r - I-kUe- D^) {ue{t,x)dx-dfit,Xx))- 
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As r — I -k fitg vanishes on the support of /itg, 



d f , . A{t) 



Ue{t,X 



- / {r - I -k Us - Ds) dfit,^ [x) = — 

^ Jr "^ Jm 

+ - / I -k{Ue- Uto) {Ue{t, X) dx - rf/iio(x)) + - / D^Ue dx, 

£ Jr ^ Jm 

with 

A{t) = / {r - I -k fis) Ue{t,x)dx+ / Ded^ito- 



Notice that 

/ -k {u^it, .) - /itj (Mj(t, x) dx - d^,s{x)) 

controls ||/ -k {us{t, .) — /ito)lli2 since P < CI. 
So we deduce since D^ > 

1 /"* f Ueit x) 

- I \\I-k{ue{s,.)- iitMl^ds< I \og-^-^dfXt^ 

^ J to JR "el'^0)-^J 

- {Ue{t,x)-Ue{to,x))dx+ ds / H,{ips{to)) dfito- 

Jr J to ^ ^ J to Jr 

(4.22) 



Step 3: Easy hounds. 
Lemma 14.11 tells that 

The total mass stays bounded in time so 

— {ui;{t,x)—Ui;{to,x))dx< {ui;{t,x)+Ui;{to,x)) dx < C. 

Jr Jr 

And 



/ log "/. ' N dfJ'to =- / {Ve{t,x) -ips{to,x))dfito 

Jr Ue{to,x) e J^ 



If 2 

< - / {ip{t, x) - ip{to, x)) dflto + - WVe - Vh^iU)- 

£ ./ro £ 
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Note that by Prop. 12. ![ /xjq is supported on {ip{s, .) = 0}. Since in addition 
we know that f < 0, 

f u,it,x) ^ 2 

/ ^°S —T—A ^/^*o <-m- ^IIl-(O). 

Consequently we deduce from (14.221) the bound 

1 /"* C 2 /"* Ais) 

- ||/*(M,(s,.)-/ii„)||i2C?s< — + -||^,-(/?||L^(n)+ / ^-^rfs. (4.23) 

Step 4- Control on A. 
The key point here is that suppyUig = {(p{to, .) = 0} is very close to {ip^ > 
—(fe} as implied by Lemma [4.111 

First of all De(s, .) is Lipschitz in x and vanishes on {ips{s, .) > — <^e}, so 
one has that 

/ / D^{s,x)d^tods <C I 5(supp/iio, {ipe{s,.)>-'^e])ds. 
Jto Jr Jto 

By the definition of /i, its support is included in {if = 0} and by Lemma H^ 
we finally deduce that 



/ / Ds{s, x) dutg ds <C T{t - to) {t - to) 

Jto Jr 



'to 

+ C [ S{Ms, .) = 0}, {^sis, .) > -^J) ds. 

Jto 

(4.24) 

This was the only additional term with respect to [9J. The other term is 
treated the same by decomposing 

/ {r - I -kij,tg)us{s,x)dxds = {r - I -k Ht^) Us{s,x)\^(^s,x)<-^,dxds 

Jto Jr JtoJR 

+ {r - I -k ^to) Ue{s,x)\^(^s,x)>-(p,dxds. 

Jto Jr 

For the first part, note again that by (12.61) . there exists R s.t. V|a;| > R, r — 
I*f^to< 0. 
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Therefore, we may simply dominate 

/ / {r - I-k fit^^) Us{s,x)I^^^s,x)<-^,dxds <C{t- s)e 

Jtn Jr 



Concerning the second part, just note that {(/^^(s,.) > — <^e} C {ip{s, .) > 
—as} with as = <^e + Hv? — v^eHloo and that r — I -k fito is non positive on 
{(p{tQ, .) = 0}. Hence 



Jto Jr 



'to 

<C [ 6{{^{s, .) > -a J, {<^(to, .) = 0}) ds. 

Jto 

By Lemma [4. 8[ this imphes that 

/ j [r -I-k /iij Ue{s, x)\^(s,x)>-g,, dx ds < T{t - to) {t - to) 

Jto JR 

+ C [ 6{{^{s, .) > -a J, {^{s, .) = 0}) ds. 

Jto 

Inequahty fl4.23p now becomes 



ft Q ^ 

/ ||/ * (m,(s, .) - Ht,)\\l2ds < ^ + 2 y^ - <^||l-(q) +C{t-s) e-^^l' 

Jto Ve 

+ CT{t- to) (t - to) + C / 5{{ip{s, .) > -«e}, {^{s, = 0}) ds 

Jto 

+ C [ 6{Ms, .) = 0}, {^s{s, .) > -(^J) ds. 

Jto 

(4.25) 

Conclusion. Eq. fl4.22p indeed gives Lemma [4.121 if one defines 
Ce 

+ C f 6{Ms,.)>-as},{^{s,.) = 0})ds 
Jo 

+ C [ 5iMs,.)=0},{ipsis,.)>-^e})ds. 
Jo 
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Thus, in order to complete the proof of Lemma 14. 12^ we only have to check 
that a{e) — )■ when e — > 0. 

It is obvious for the first three terms. For the second note that 

C 5{{(p{s,.)>-ae},{(pis,.) = 0})ds — ^0 as £ -^ 0, 
Jo 

since by dominated convergence it is enough that for any s 

<5({<^(s,.)>-aJ,{^(s,.) = 0})^0. 

For the last term, we use Lemma |4]TT] to show that it converges to as e — > 0. 
Indeed for any 6 one has 

/ 6{{^{s,.) = 0}, {^,{s,.)>-ip,})ds<\J,\+T6, 
Jo 

with I Je| — !• 0. □ 

4.7 Concluding Remarks: Finite speed of propagation 
of the support 

It is now easy to obtain additional properties for the limit ip, for instance the 
finite speed of propagation of the support. We just sketch quickly here how 
such properties can be proved. 

First by the definition of /i, one has that 

r-J^/i({(^(t,.)})<0, on{(^(t,.) = 0}. 

One then deduces immediately that for any x 

r{x)-I^f^{{ip{t,.)})-Dit,x)<0, 

since D = m.m{K d{x, {ip(t, .) = 0}), Dq). That means by Eq. (14.211) 

dt^p < H{d^(p). 

Just by the following the characteristics of the Hamilton- Jacobi equation, 
one obtains that the set {(p = 0} can only propagate at finite speed V. 

Finally let us mention that some time regularity is also known on /i. 
For instance the set {cp = 0} is continuous in time by Lemma 14.21 By the 
stability result in p3], one deduces that fi{{ip{t, .)}) is continuous in time. 
As a matter of fact the finite speed of propagation of {ip = 0} would even 
imply Holder continuity of order 1/2 for /i. 
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